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ABSTRACT 

A series of optical and one near-infrared nebular spectra covering the first year of the Type la 
supernova SN 201 Ife are presented and modelled. The density profile that proved best for the 
early optical/ultraviolet spectra, “p-1 Ife”, was extended to lower velocities to include the re¬ 
gions that emit at nebular epochs. Model p-l Ife is intermediate between the fast deflagration 
model W7 and a low-energy delayed-detonation. Good fits to the nebular spectra are obtained 
if the innermost ejecta are dominated by neutron-rich, stable Fe-group species, which con¬ 
tribute to cooling but not to heating. The correct thermal balance can thus be reached for the 
strongest [Fe ll] and [Fe III] lines to be reproduced with the observed ratio. The ^®Ni mass thus 
obtained is 0.47 ± 0.05 Mq. The bulk of ®®Ni has an outermost velocity of ~ 8500km s“^. 
The mass of stable iron is 0.23 ± 0.03 Mq. Stable Ni has low abundance, ^ 10“^ Mq. This is 
sufficient to reproduce an observed emission line near 7400 A. A sub-Chandrasekhar explo¬ 
sion model with mass 1.02 Mq and no central stable Fe does not reproduce the observed line 
ratios. A mock model where neutron-rich Fe-group species are located above ®®Ni following 
recent suggestions is also shown to yield spectra that are less compatible with the observa¬ 
tions. The densities and abundances in the inner layers obtained from the nebular analysis, 
combined with those of the outer layers previously obtained, are used to compute a synthetic 
bolometric light curve, which compares favourably with the light curve of SN 201 Ife. 

Key words: supernovae; general - supernovae: individual (SN 201 Ife) - techniques: spec¬ 
troscopic - radiative transfer 


1 INTRODUCTION 

Given the importance of Type la supernovae (SNe la), one fun¬ 
damental question concerns their origin. There is strong evidence 
that they are the explosion of compact stars, most likely carbon- 
oxygen (CO) white dwarfs (WDs). Such stars can undergo ther- 
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monuclear explosive burning which results in the production of 
radioactive ®®Ni, whose decay powers the SN luminous output, 
and other elements along the a-chain. Apart from ®®Ni, the ele¬ 
ment which is produced most abundantly is silicon, whose pres¬ 
ence characterises the spectra of SNe la at early times jNugent et alj 
l201lLlFolatelli et al.ll2012l:|Parrent et ahkoia We kn ow that ®®Ni 
is produced in variable amounts jMazzali et alj|20(y^. and that the 
relation between ^^Ni, lu minosity, and opacitv lpinto & EastmanI 
120001 : [Mazzali et alj200lh shapes the light curves of SNe la, which 
are to first order a one-parameter family where b roader light curves 
belong to more luminous SNe jPhilliDslll993h . The explosion is 
thought to start when the WD reaches a central temperature of 
~ 10® K, which happens when its mass approaches the Chan¬ 
drasekhar limit, or it may be triggered by accretion or merger events 
in lower-mass WDs. 

What we do not know for sure, despite much effort, is how 
the CO-WD reaches such a large mass, and how it then explodes. 
We need to establish this so that we can then tackle the problem 
of how the explosion can produce SNe with different ®®Ni content 
and hence different luminosity. Popular models for pushing the WD 
mass close to the Chandrasekhar limit involve transfer of mass from 
a nondegenerate binary companion (main sequence or red giant), 
the single-degenerate (SD) scenario, or the merger of two WDs, 
the double-degenerate (DD) scenario. Both scenarios have pros and 
cons. 


The SD scenario ensures repeatability, a constant exploding 
mass, and a spread of ages, but it is supposed to involve hydro¬ 
gen accretion, which should lead to the detection of H, either in 
the outer layers as a leftover of the accreting material or at low 
velocities as i s expected of H stripped from the com panion by 
the SN ejecta ([Marietta et al.ll200ol : IPakmor et al.ll2oTIh. However, 
only r arely has H been directly detected in SNe la dSilverman et all 
l2013h . Interes tingly, in these cases much H is s een, such as 
in SN2002ic dHamuv et al.ll2003l: jPeng et al.ll2004 or PTFllkx 
dPildav et al.ll2012l) . Also, the identification of a surviving com¬ 
panio n is not certain llRuiz-Lapuente et '^120041 : iKerzendorf et al.l 

I2OO9I1 . 

The DD scenario may occur more frequently. It is, however, 
unclear how many merger events actually exceed a combined mass 
of about the Chandrasekhar m ass. The prop osed scnearios, violent 
mergers i PakmoreLaL 201T) o r colli sions dRosswog et al.l 1^091 : 
lThompsonll201 ll : iKushnir et al.ll2013l) . produce diverse outcomes 
depending on the initial conditions. The question then is whether 
they look like a typical SNe la, especially in view of the fact that 
they involve large asymmetries because ignition happens at the 
point of contact rather than at the highest density. 

In a further twist, the possibility that WDs below the Chan¬ 
drasekhar limit also explode if they accrete He-rich matter from a 
companion has also been proposed dLivne & ArnetlllT995h . This, 
the sub-Chandra scenario, may again lead to a different set of 
obser ved properties (e.g., [Nugent et alJll99% IWooslev & KasenI 

IMIh . 


Since all of these scenarios are based on solid physical ideas, 
it is perfectly possible that they are all actually verified in nature. 
Is this indeed the case? Does this add to the uncertainty in our use 
of SNe la as distance indicators? While the latter is a question that 
can also be addressed by using large samples of SNe, the former is 
one that requires the accurate study of a few, well-observed SNe. 

The biggest differences among explosion models are expected 
to be found either at high velocities (presence of circumstellar mat¬ 
ter, amount of material at high velocity as an indication of kinetic 
energy, metal content of the progenitor, unbumed material from the 


progenitor or the companion), or at the lowest ones (central densi¬ 
ties, composition, material dragged from companion upon impact). 
Different epochs in a supernova’s evolution yield different infor¬ 
mation. Early-time data can be used to infer the size of the progen¬ 
itor, possibly that of the companion (in the SD scenario), and the 
properties of the outer layers of the exploding star, which may not 
be affected by thermonuclear burning and may thus carry informa¬ 
tion about the accreting material. Observations over the first month 
or so tell us about general properties: the light-curve shape (mass, 
opacity), ejecta velocity (kinetic energy), and composition of the 
outer layers (extent of ®®Ni, intermediate-mass elements (IMEs)). 
Correlations have been found among these quantities, appearing to 
suggest some regularity. 


However, if information about the central regions is required, 
one has to wait for the nebular phase. At late times 1 yr), opac¬ 
ities are low and it is possible to explore the lowest-velocity mate¬ 
rial. This corresponds to the highest densities, and hence it should 
bear the trace of the ignition and the outcome of the explosion. 


SN 201 He ( [Nugent etal.l201lh is the perfect candidate for ac¬ 
curate studies. When it was discovered it was the nearest SNIa 
in almost 40 yr. It h as been observed e xtensively at frequen - 


from the radio ([Horesh et al.[ [20121 : Chomiuk et alJ 2012h , 


to the far-infrared (IR) with Herschel ([Jo hansson et alj 2013[) . to 


iMatheson et alj|2012h. to the optical jNusent et alJ 

201 ll: Li etal. 

201 ll: Richmond & Smith[[2012l:[Munari et al.[[2013l: 

Shannee et al. 


2 OI 3 I : Pereira et al. 2 OI 3 I: Tsvetkov et alj 2013 ), to the ultra vio- 


let (UV) with Swift ([Nugent et al 


with the HST ( [Foley & KirshneJ 


[2011[:Ib 


rqwnet 




2013[: [Mazzali et al 


2012h and 


[2014. to 


the X-rays w ith Swift 1 Hore sh et alj 2 OI 2 I : Margutti et ^ 2012ll 


and Chandra ([Margutti et al.[[2012[) . and finally to the gamma rays 
([isem et al.[[2013). 


SN2011fe was a normal SNIa hosted in a spiral galaxy, 
MlOl. The early discovery of the SN by the PTF collabora¬ 
tion allo wed a good estim ate to be m ade of the time of ex- 
plosion ([ Nugent e t alj [201 iL revised by [Piro & Nakail [20131 and 
[Mazzali et al.[ 2014h . and also made it possible to estimate the 
size o f the exploding star to be R^mg 0.02 Rq ([Bloom et alJ 
[2012h . The early-time light curve did not reveal the signature ex¬ 
pected from a shoc ked companion, con t rary to predi ctions if this 
were a sizable star ( [Nugent et alj[20lll : [KasenIl2oTc ). Early-time 
spect r a exhibit traces of c arbon in the outer layers 1 Parrent et aH 
[20 12I : [Mazzali et al.[[2014[) and are consistent with the low metal- 
licity of the host galaxy ([Mazzali et alj[2014h . The spectra have 
been claimed to be con sistent with both a SD and a DD sce¬ 
nario ([Rbnke et al.[[2012h . In the outer layers, a density interme¬ 
diate between the fast deflagration model W7 and a low-energy 
delayed detonatio n was found to be t he best solution for the op- 
ticalAJV spectra ([Mazzali et alj [2oI4 . Merger models and sub- 
Chandrasekhar-mass models seem to have lower densities in the 
inner layers, and a nebular study can help determining how well 
they might reproduce the properties of SNe la. 


In the following, the nebular-phase data for SN 201 He are pre¬ 
sented (Section 2). We show models of the optical spectra in Sec¬ 
tion 3 using the density profile derived from the early-time spectra. 
In Section 4, we extend the models to the near-infrared (NIR) and 
discuss alternative models. A discussion (Section 5) and conclusion 
(Section 6) follow. 
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2 NEBULAR SPECTRA OF SN 2011FE 


Table 1. Observing log for the nebular spectra of SN 201 Ife 


A series of 10 nebular spectra of SN 201 Ife was taken over the 
period March to August 2012. An observing log can be found in 
TableE 

Four spectra were taken with the William Herschel Telescope 
(WHT) using the Intermediate dispersion Spectrograph and Imag¬ 
ing System (ISIS). ISIS is a dual-arm spectrograph; we used the 
R158R grating and GG495 order hlocker (red arm), as well as the 
R300B grating (blue arm) and the 5300 A dichroic, giving a wave¬ 
length coverage of ~ 3500-10,000 A. Four spectra were taken us¬ 
ing the Lick Observat ory 3-m Shane telesco pe and the Kast dual¬ 
channel spectrograph l lMiller & Stonelfr993h . The 600/4310 grism 
was used in the blue arm and the 300/7500 grating in the red arm, 
giving a wavelength coverage of 3450-10,200 A. 

All of t he Lick spectra were reduced using standard tech¬ 
niques (e.g., ISilverman et al.ll2012h . Routine CCD processing and 
spectrum extraction were completed with iraiQ We obtained the 
wavelength scale from low-order polynomial fits to calibration- 
lamp spectra. Small wavelength shifts were then applied to the 
data after cross-correlating a template sky to the night-sky lines 
that were extracted along with the SN. Using our own IDL rou¬ 
tines, we fit spectrophotometric standard-star spectra to the data 
in order to flux calibrate our spectra and to remove telluric lines 
dWade & Homelll988l;&heson et aljioooh . 

Spectra of SN 201 Ife were obtained with the Large Binocu- 
lar Telescope (LB T) and Multi-Ohject Dual Spectrograph (MODS; 
jPogge et^l201 2ll on three nights in 2012: March 24, April 27, 
and June 12 (UT dates are used throughout this pa per). Details of 
the da ta acquisition and reduction can he found in IShappee et al.l 
( l2013h . 

NIR spectra were also o btained with the LBT LUCIFER spec¬ 
trograph dSeifert et alj|2003l) on 2012 May 1, close to the date of 
one of the MODS optical spectra. We used the 210 line mm~^ grat¬ 
ing and 1.0"-wide slit, resulting in a resolution of 4000 in the J 
hand. Data were taken at three grating tilts to obtain spectra in the 
.7, 77, and K bands. Consecutive exposures were dithered along 
the slit and used to subtract the sky background. Images at each slit 
position were then combined and extracted using the IRAF “twod- 
spec” package, and the two spectra were averaged to create the final 
spectrum. The spectra were wavelength calibrated using night-sky 
emission lines, and telluric absorption features were removed using 
a spectrum of the AO V star HIP67848 taken immediately after the 
SN spectra. The NIR magnitudes of HIP67848 were then used to 
flux calibrate the spectra of SN201 Ife. 

An acquisition image of SN 201 Ife was obtained in the J band 
with LUCIFER. Using two 2MASS catalog stars in the field we 
estimate the brightness of the SN to be J = 17.20 ± 0.06 mag. 

The spectra were reduced using custom-written pipelines 
based on standard procedures in IRAE and IDL, broadly follow¬ 
ing the reduction procedures outlined bv lEllis et al.l ( 1200 8h . includ¬ 
ing flux calibration and telluric- feature removal. A ll spectra were 
observed at the parallactic angle dFilippenkdl 1 982h in photometric 
conditions. “Error” spectra are derived from a knowledge of the 
CCD properties and Poisson statistics, and are tracked through¬ 
out the reduction procedure. The spectra are also rebinned (in a 


^ IRAF is distributed by the National Optical Astronomy Obsei'vatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation (NSF). 


UT Date 

Telescope 

Instrument 

Wavelength 

Coverage 

2012 Mar 2 

WHT 

ISIS 

3500-10,000 A 

2012 Mar 24 

LBT 

MODS 

3300-10,000 A 

2012 Apr 2 

Lick 

Kast 

3450-10,300 A 

2012 Apr 23 

Lick 

Kast 

3450-10,200 A 

2012 Apr 27 

LBT 

MODS 

3300-10,000 A 

2012 May 1 

LBT 

Lucifer 

1.17-1.31 pm 

2012 May 1 

LBT 

Lucifer 

1.55-1.74 pm 

2012 May 1 

LBT 

Lucifer 

2.05-2.37 pm 

2012 May 26 

WHT 

ISIS 

3500-9500 A 

2012 Jun 12 

LBT 

MODS 

3300-10,000 A 

2012 Jun 25 

WHT 

ISIS 

3500-10,600 A 

2012 Jul 17 

Lick 

Kast 

3450-10,300 A 

2012 Aug 21 

WHT 

ISIS 

3500-10,000 A 

2012 Aug 23 

Lick 

Kast 

3450-10,200 A 


weighted way) to a common dispersion, and the two arms of the 
spectrographs joined. 

The nebular spectra of SN 201 Ife have been checked ag ainst 
the late-time broad-band photometry of iTsvetkov et alj (120131) . The 
broad-band magnitudes at the epochs of the nebular spectra were 
derived, for each band, from a linear fit of the available photometric 
points on the exponentially declining tail of the light curve. The flux 
calibration of the spectra was then checked against the photome¬ 
try using the IRAF task STSDAS.HST.CALIB.SYNPHOT.CALPHOT. 
Next, we locked the absolute flux scale of each spectrum, normal¬ 
ising the synthesised BVRI magnitudes from the spectrum to the 
interpolated broad-band magnitudes. The 2012 Aug. 21 and 2012 
Aug. 23 spectra were combined into a single spectrum to improve 
the signal-to-noise ratio. 


3 MODELLING 


We modelled the nebular spectra of SN201 Ife using our well- 
tested SN n ebular code. Th e code is based on approximations de¬ 
veloped by I Axelrod (Il98(lh. and it has been developed and dis - 
cussed in various papers ( Mazzali et al.ll200i] ; lMazzali et al.l2007h . 
The SN ejecta are divided into a number of spherical shells, whose 
density reflects that of an explosion model and whose composition 
can be varied arbitrarily. The ejecta size and density are rescaled 
to the epoch of the desired model assuming homologous expan¬ 
sion. ®®Ni decay provides the heating. In the first step of the pro¬ 
cedure, 7 -rays and positrons are emitted according to the distribu¬ 
tion of ^ ^Ni. Their propagation is followed with a Monte Carlo 
scheme dCappellaro et alj 1997h . While 7 -rays can travel a long 
way before being absorbed in the late-time low-density SN ejecta, 
positrons mostly deposit their energy locally. Deposition of the en¬ 
ergy carried by 7 -rays and positrons determines the heating rate. 
Under nebular conditions, this collisional heating is balanced by 
cooling via line emission. The ionisation and excitation of relevant 
species is computed assuming non-local thermodynamic equilib¬ 
rium (NLTE), balancing heating and cooling in an iterative scheme. 
When convergence is reached, line emissivity is transformed into 
an emission-line spectrum. An emitted spectrum is computed for 
every shell, and the final emerging spectrum is the sum of the con¬ 
tributions from all shells. 

A successful model should be able to reproduce the evolution 
of the spectrum with time, reflecting the decreasing heating rate and 

































4 P. A. Mazzali et al. 


Table 2. Nebular models 


Date 

Epoch 
[if days] 

M(CO) 

Mq 

M(IME) 

Mq 

M(5®Ni) 

Mq 

M(54Fe) 

Mq 

M(®®Ni) 

Mq 

M(NSE) 

Mq 

i;(l-zone) 
km s~^ 

2012 Mar 2 

192 

0.24 

0.42 

0.41 

0.29 

0.007 

0.71 

8000 

2012 Apr 2 

223 

0.24 

0.42 

0.44 

0.26 

0.007 

0.71 

7900 

2012 Apr 23 

244 

0.24 

0.42 

0.45 

0.25 

0.006 

0.70 

7900 

2012 Apr 27 

248 

0.24 

0.42 

0.46 

0.24 

0.004 

0.70 

7800 

2012 May I 

252 

0.24 

0.41 

0.48 

0.22 

0.004 

0.69 

7800 

2012 May 26 

277 

0.24 

0.41 

0.46 

0.22 

0.007 

0.70 

7900 

2012 Jun 25 

307 

0.24 

0.41 

0.48 

0.21 

0.007 

0.70 

7600 

2012 Jul 17 

329 

0.24 

0.39 

0.51 

0.22 

0.007 

0.73 

7500 

2012 Aug 22 

364 

0.24 

0.38 

0.54 

0.18 

0.008 

0.72 

7700 


deposition. iMazzali et al.l ( |2014|) showed that the early-time spec¬ 
tra of SN 201 Ife cannot be satisfactorily reproduced with available 
SN la explosion models, and proposed a density distribution char¬ 
acterised by a small but nonzero mass at high velocities, above 
~23,000kms“^ (see Fig. 0. This model was shown to be able 
to reproduce both the optical and UV spectra simultaneously. The 
early-time study could not explore the deepest parts of the ejecta, 
which are optically thick at early epochs. Flere we extend the re¬ 
sults of iMazzali et al.l ( l2014h and explore the deepest regions of the 
ejecta. 

iMazzali et alj j2014h concluded that in SN 201 Ife some ~ 
0.4 Mq of ®®Ni are located at u > 4500 km s“^. Estimates 
of the t otal ^^Ni mass for SN 201 Ife suggest a value close to 
0.5 Mq jPereira et al. I l2013h . If correct, this implies that only a 
small fraction of all ®®Ni is located at low velocities. If the total 
ejected mass of SN 201 Ife was close or equal to the Chandrasekhar 
mass, using the density structure p-llfe leaves 0.25 Mq below 
4500 km s “ ^. Additionally, because of limitations of our early-time 
code caused by the lower boundary black-body assumption (e.g., 
IMazzali et al.lll99^ . abundances in layers where ®®Ni is dominant 
may not be fully reliable. In SN 201 Ife this affects the region be¬ 
low ~ 7000 km s~^. In an y case, we adopted the abundances de¬ 
termined by IMazzali et al.l ( |20I4|) as our initial distribution, and 
adapted those in the region below 4500 km s“^ in order to fit the 
spectra. The density structure of model p-llfe, as well as those of 
other models used in this paper, are shown in Figure [T] 

Assuming that the spectra are properly calibrated in flux (there 
may be some uncertainty when spectra were taken at different tele¬ 
scopes, see below), the overall flux level should reflect the to¬ 
tal mass of ®®Ni. An important remark here is that even in the 
case of SN 201 Ife the uncertainty in flux calibration can account 
for an uncertainty in the luminosity level of up to 10 per 
cent. In the modelling, we adopted for SN 20 1 Ife a distance of 
6.4 M pc (distance modulus p = 29.04 mag: IShappee & StanekI 
I2OI iL find p = 29.040.05 (random) 0.18(systematic) mag), and 
reddenings E{B — V) = 0.014 mag in the host galaxy and 
E(B — V) = 0 .009 mag in the Milky Way iSchlegel et aljf 199i : 
[Patat et alj(2013l) . 

We find that the nebular spectra of SN 201 Ife require on av¬ 
erage a ®®Ni mass of ~ 0.47 Mq. This is con sistent with previou s 
estimates based on the peak of the light curve dPereira et al.ll2013h . 
Our set of nebular models is shown in Figures and and Ta¬ 
ble 1 summarises the main results. There seems to be a tendency 
for the ®®Ni mass to increase for later-time spectra. This may be 
caused by inaccuracies in our description of positron escape, sinc e 
little is known about positron deposition dCappellaro et al.l[l997h . 
but in test models with full positron trapping we are getting similar 



Figure 1. Density structures used in this paper: W7 (blue), delayed- 
detonation model WDDl (light blue, solid), sub-Chandra model CDT3 
(grey), p-llfe (red), and a sub-Chandra p-llfe model with M = 1.2 Mq 
(light red). 

results. Another possibility is that too much IR flux is produced by 
our models for the later epochs. There may be indications of this in 
the Fe-dominated emission features in the red, which are increas¬ 
ingly overestimated by the model spectra, while the emission at 
9400-9800 A is not reproduced in our models. These are high-lying 
Fe II transitions. One possibility is that the collision strengths of 
these and other weak lines are not well determined. Overall, how¬ 
ever, given the poor fitting of some features, possible uncertainties 
in flux calibration, and the neglect of radiative-transfer effects, we 
regard our ®®Ni mass estimate to be within the uncertainties. 

Another region that is poorly reproduced is the emission near 
4250 AA. Several moderately strong [Fe II] lines are present there, 
including lines at 4244 AA, 4277 AA, and 4416 AA, but their com¬ 
bined strength is only about half of what is observed. This may 
be due to problems with collision strengths, A-values, or it may 
indicate the presence of some other species which is not consid¬ 
ered in our models. The same is true for the narrow emission near 
5000 AA. 

The element distribution determined from the early-time spec¬ 
tra is compatible with the late-time models. These are really only 
sensitive out to u « 8000 km s“^, which is the width of the emis¬ 
sion lines. Fitting the Fe emission lines, we find a roughly constant 
width, perhaps with a slight tendency to decrease with time. The 
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rest wavelength [A] 


rest wavelength [A] 


Figure 2. Synthetic spectra based on the p-1 Ife model (blue) compared with observed spectra (black): earlier spectra. 


characteristic width at epochs of ~ lyr is ~ 7600 km s“^. Based 
on a relation betwe en emission-line width and light-curve shape 
dMazzali et aljl99^ . we expect for SN 201 Ife a full width at half¬ 
maximum intensity of the strongest Fe line blend of ~ 12,500 km 
s“^, and this is indeed what we measurj^ This corresponds to a 
normal SNIa, but appears to be somewhat narrow for the ®®Ni 
mass and the light-cur ve shape (Ami 5 {B) = 1.07 ± 0.06mag; 
[McClelland et al.ll2013h . 


^ The width of the emission line reflects both velocit y broadening and th e 
blending of several lines that contribute to the feature iMazzali et al .Il998h . 


Our total ®®Ni mass is only slightly larger than the mass de¬ 
termined from our early-time analysis at velocities v > 4500 km 
s"*^. Only an additional ~ 0.1 M© of ®®Ni is found to be located at 
V < 4500 km s”*^. As we mentioned, model p-llfe has ~ 0.25 M© 
of material below 4500km s“^. This means that ~ 0.15 Mq of 
material at low velocity is not ®®Ni, but rather stable nuclear statis¬ 
tical equilibrium (NSE) species. Two major species are predicted 
to be synthesised in n-rich burning: ®^Fe and ®®Ni, both of which 
are stable isotopes. We have a direct handle on the abundance of 
stable Ni from the emission line near 7400 A, which is normally 
identified as [Ni II] AA7380, 7410. Indeed, the small emission seen 
in SN 201 Ife near 7400 A can only be reproduced if some Ni is 
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Figure 3. Synthetic spectra based on the p-llfe model (blue) compared with observed spectra (black): later spectra. 


present at the lowest velocities. Given the advanced epoch, this 
must be stable Ni, presumably ®®Ni. This is produced from ®‘*Fe 
in the so-called alpha-rich freezeout phase via ®^Fe {a, 7 ) ®®Ni 
during the early expansion. However, the broad emission in that re¬ 
gion is dominated by other lines, [Fe II] AA7388, 7452 and [Ca II] 
AA7291, 7324. Only a very small mass of ®®Ni is required to re¬ 
produce the small emission that corresponds in wavelength to the 
[Ni II] lines. On average, a value of ~ 0.008 M© of ®®Ni, confined 
below V « 4000 km s“^ and concentrated at the lowest velocities, 
is sufficient. Even at the lowest velocities, the abundance of ®®Ni 
never exceeds 10%. This means that most of the mass at low veloc¬ 
ities is composed of Fe. Again, this Fe cannot all come from ®®Ni 


decay, because if it did the SN luminosity would be too high. We 
presume that this is mostly ^'^Fe. 


Although we cannot distinguish ®^Fe from directly synthe¬ 
sised ®®Fe, the latte r is supposed to be p roduced at low Ye, in¬ 
side the ®®Ni region jlwamoto et al.|[l99^ . which is not where we 
find it t o be located in our modelling. This is a common feature in 


SNela dStehle et alj|2005l: iMazzali et al]|2008l: iTanaka et al.ll201ll: 
ISasdelli et al.l 201 ^ . and it may indicate a higher Ye than predicted 
by models such as W7. The presence of stable Fe helps keeping 
the ionisation degree sufficiently low that the correct relative con¬ 
tribution of [Fe II] and [Fe III] lines is achieved in our models. 
We showed in other work that if this central density is low, ion- 
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Figure 4. Abundanc e distribution in the ejecta of SN201 Ife. The elements that are shown are those that can be inferred from spectroscopic modelling of 
iMazzali et al.l I2014II . The top panel is linear with respect to mass, the bottom panel with respect to velocity. 


isation becomes too high and [Fe lllj lines do minate, leading to 
a strong 4800 A emission (IMazzali et alJIlOllh . If the density is 
even lower, only [Fe 111 ] lines are produced and the emission-line 
spectrum changes nature. This seems to be a f eature of the fastest- 
declining SNela jMazzali & Hachingeill20l3l . The total mass of 
stable Fe (o btained by adding m aterial already diagnosed at high 
velocities by IMazzali et alj[20l4 which is necessary to reproduce 
Fe lines at early times, when ®®Ni has not yet decayed to ®®Fe, 
and stable Fe at low velocities as determine d here) is on average 
~ 0 .23 Mfn. This is in line with the results of iMazzali et al.l l l2007ll 
and iKasen & WooslevI ( 120071) . 


A plot of our abundance distribution is shown in Figure|4] This 
compares favourably with one-dimensional (ID) Chandrasekhar- 
mass models based on the SD scenario, where ignition occurs near 
the centre of the WD at high density, and a sufficiently high tem¬ 
perature is achieved that n-rich isotopes can be synthesised. 


3.1 The Near-Infrared 

We mentioned that as time goes on, the flux is expected to shift to 
the infrared. We do have one epoch of NIR spectroscopy, obtained 
with the LET on 1 May 2012. This spectrum can be combined with 
the optical spectrum obtained on 27 April 2012. In Figure we 
show our model for that epoch, now extending to the NIR. The 
model reproduces the J- and JF-band spectrum reasonably well, 
while the flux in the H band is too low. This is in part because [Si I] 
1.6 /im is overestimated, the result of ®®Ni being in close contact 
with Si, but we also seem to somewhat overpredict the Fe emission 
in that region. 

It is possible that a sharper velocity sep aration between ^^N i 
and IMEs may be required than the results of iMazzali et alj | |2014|) 
suggest. We tested this in a model, but the resulting spectrum still 
showed too much flux at 1.6/rm. Using a three-dimensional (3D) 
approach may be useful, concentrating ®®Ni in blobs or clumps 
from which 7 -rays escape less successfully. On the other hand, 
while this may help, 7 -rays escape quite easily at the low densities 
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Figure 5. Synthetic spectrum based on the p-llfe model (blue), now also showing the IR. The observed IR spectrum (1 May 2012) is given along with the 
optical part (27 Apr 2011). The strongest lines in the model, as well as lines that match observed features in wavelength but possibly not strength, are labelled. 


that in general characterise the nebular regime, so the overdensity 
would probably have to be quite large. 

The other possibility is that the density in the inner region may 
still be somewhat too high. In the case of SN 2003hv, a sharp reduc¬ 
tion of the NIR flux was obtained when the central densities were 
decreased, shifting the flux from [Fe II] to [Fe III] lines. The ratio of 
these two ions in the optical is correctly reproduced for SN 201 Ife, 
however, so we do not expect that the density can be reduced by 
a large amount. In any case, we can test this possibility using sub- 
Chandrasekhar-mass models. 


3.2 sub-Chandrasekhar-mass models 


In order to verify the need for a core of stable NSE elements, we 
tested a sub-Chandrasekhar-mass model. In such models, explo¬ 
sive burning begins at the outside, typically following the accre¬ 
tion of He-rich material from a companion. Most models suggest 
that a secondary detonation can be triggered within the WD, which 
can be incinerat ed even if its mass is s i gnificantly below the Chan¬ 
drasekhar limit jLivne & Amett|[l9^ : iFink et alJbOlOh . However, 
because the density never reaches the high values typical of near- 
Chandrasekhar-mass WDs, n-rich nucleosynthesis is not expected 
to take place. Models are characterised by low central densities. 


We used the density and abundance distribution of a 
ID sub-Chandrasekhar -mass model from the set presented by 
IShigevama et al.l l ll992t) . We selected model CDT3, which yields 
a ®®Ni mass of 0.50 M©, in agreement with the estimated value. 
The total mass of the model is IM©. The mass of stable 
Fe-group isotopes, ®^Fe and ®®Ni, is very small, 0.003 M© 
each. Intermediate-mass elements are present with significant abun¬ 
dances: Si, 0.21 M© and S, ~ 0.12 M©. Some unbumed oxygen 
is also present, ~ 0.07 M©. The density structure of this model is 
shown in Figure[T] 

The corresponding synthetic spectrum computed for day 248 
is shown in Figure]^ It is characterised by high ionisation, as shown 
by the ratio of [Fe III] and [Fe II] lines: the [Fe III] -dominated 
emission near 4800 A is much stronger than the [Fe II] -dominated 
emission near 5200 A. This is reminisce nt of SN 2003hv, and 
is the result of th e low central density dLeloudas et al.l 1 20091 : 
iMazzali et alj|20li]) . Elsewhere, the small Si mass, its separation 
from ®®Ni, and the lack of Fe II means that the synthetic line at 
~ 1.6 irm is now weaker than the observed one. However, other 
synthetic [Fe II] lines, such as those near 8800 A, are still too 
strong, suggesting a problem with collision strengths or, less likely, 
^-values. On the other hand, the [Ni II] line is no longer repro¬ 
duced. Also, very noticeably, the Ca abundance is much too high. 
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Figure 6. Sub-Chandra models: CDT3 (blue) and a model with M = 1.26 Mq (red). The observed spectrum is shown in black. 


as shown by the strength of the emission near 7200 A. This demon¬ 
strates that a model which is significantly sub-Chandrasekhar is not 
appropriate for SN 201 Ife. It is much more difficult to establish 
whether a model which is just below the Chandrasekhar mass can 
also be ruled out. This mostly depends on the exact nucleosynthe¬ 
sis. It also confirms that stable Fe is required in order to reproduce 
the spectra. 

As an alternative, we constructed a sub-Chandrasekhar model 
with a larger mass. In this test model the total mass is 1.26 Mq, 
the stable Fe content is small but not zero (0.12 M©), and the ®®Ni 
mass is 0.43 Mq, consistent with our results above. The synthetic 
spectrum (Figure shows again the problem of high ionisation, 
although to a lesser extent. The best results is obtained when the 
mass approaches Mch and the stable Fe content is close to 0.2 Mq. 

3.3 Where is the stable Fe? 

Having established that a high central density is needed, and that 
the presence of a significant amount of ®'*Fe and stable ®®Fe pro¬ 
duced at burning l llwamoto et al.l[l99^ appears to be necessary, 
we can address the question of its location. Recent 3D delayed- 
detonation models suggest that n-rich NSE species produced in the 
deflagration phase rise because of their buoyancy and end up ex¬ 
panding at higher velocities than ®®Ni when the simulation is per¬ 


formed in 3D. This leads to an i nverse structure, where ®®Ni is lo¬ 
cated inside n-rich NSE species dSeitenzahl et alj|20 1311^1 In order 
to assess the effect of this, we simulated spectrum formation us¬ 
ing the p-1 Ife density structure but moving stable Ee and Ni above 
®®Ni in velocity. This process requires a moderate rescaling of the 
®®Ni mass to obtain the same luminosity, because when ®®Ni is lo¬ 
cated deeper in the ejecta 7 -rays deposit more efficiently. The flux 
level of the nebular spectrum can be reproduced for a ®®Ni mass 
of only 0.34 Mq. The mass of stable Fe becomes 0.43 Mq, and 
®®Ni slightly increases to 0.01 Mq. The resulting model is shown 
in Figure[7] 

Although the model has reasonable features, three areas cause 
concern. First, all Fe lines are narrow. This is because there is too 
little ®®Ni at high velocity, and the Fe that is located there is not 
efficiently heated by the 7 -rays, which are emitted very deep in 
the ejecta. Increasing the mass of ®®Ni would help, but the flux 
would be too high. Additional mixing may fix this problem, but 
it may take us back toward a situation where ®^Fe is below ®®Ni. 
Second, the ratio of Fe lines now favours Fe II. This is not sur¬ 
prising, since radiation comes predominantly from the densest, in- 


® This did not see m to have been noticed in earlier 3D simulations 
iOamezo et alj200^ . 
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Figure 7. Synthetic spectrum for the Ni-Fe inversion model (blue). The observed spectrum is shown in black. 


nermost regions, where recombination is most effective. The ratio 
of the two strongest emission lines, whose atomic data are proba¬ 
bly best known, is a driving parameter to establish the success of a 
model. Finally, the ®®Ni mass is in rather strong disagreement with 
the value derived from the peak of the light curve. Combining the 
smaller ®®Ni mass with the fact that the diffusion time of the optical 
photons can only increase if they are produced deeper in the ejecta, 
as would be the case if ®®Ni were located at the lowest velocities, 
this ad-hoc model is not expected to reproduce the observed light 
curve. 

In conclusion, the nebular spectrum is not as well reproduced 
by this model as it is when ®®Ni is located outside stable Fe. Also, 
whether such a small ®®Ni mass is able to reproduce the peak of the 
light curve is highly questionable. We tend not to favour this mock 
model. Tests should be performed with the real model, which is 
obviously technically superior to the ID calculation. 


4 LIGHT-CURVE MODEL 

The description of the innermost part of the ejecta deduced through 
models of the late-time spectra can be combined with that of 
the outer layers obtained through the analysis of the early-time 
spectra. Figure |4] gives the abundance distribution, while the den¬ 


sity profile was shown in Figure [T] Although the nebular results 
seem to confirm the distribution in density that was derived from 
the early-time data, the early- and late-time analyses are some¬ 
what disconnected. A test for the overall model, in the spirit of 


s pmt ot 

Ill 2008 l: 


abundance tomography l Stehle et ^ 20051 : iMazzali et alj _ 

iTanaka et alj|201 ll ; ISasdelli et alj|2014h . can be provided by the 
independent study of the light curve. The density and abundance 
stratification can be used to compute the emission and deposi¬ 
tion of 7 -rays and positrons from ®®Ni and ®®Co decay, the prop¬ 
agation of optical photons, and the final emission of radiation. 
We em ploy a Monte Ca r lo cod e developed by P.A.M. and first 
used by I Cappella ro et al.l l ll997h . The opacity description follows 
IMazzali et alj l 2001), with different weig hts assigned to different 
element groups (see lKhokhlov et al.l ll993 ) in order to descr ibe the 
dominance of line opacity in a SN la jPauldrach et all 19961) . 


The resulting bolometric light curve is shown in Figure [8] 
There we also give a pseudo-bolometric light curve computed 
based on the available UV through IR photometry. It is very sim¬ 
ilar to the light curve published by [Pereira et al.l d2013h in the re¬ 
gion of overlap but it extends to later epochs. The epochs of the 
observed light curve were of fset to match the explosion date de¬ 
rived by IMazzali etaP j2014|) — 22 Aug. 2011, MJD55795.3 — 
implying a rise time of 18.3 days. The agreement between the syn¬ 
thetic light curve and the observed one is quite good, both in flux 
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Figure 8. Synthetic bolometric light curves based on the density of p- 
llfe a nd the abundance stratification derived here and by iMazzali et alj 
(2Q14t). The UV th rough IR pseudo-bolometric light curve computed by 
IPereira et al ] i2on is shown as black circles and ours as red circles. The 
uncertainties in the flux are as quoted in that paper. An equivalent bolo¬ 
metric light curve computed by S.B. is shown as red circles. The epoch of 
both bolometri c light curves has be en shifted to match a risetime of 18.3 
days, following lMazzali et al.l feOldh . The lines show two light curves com¬ 
puted with our density and abundance results: = 7 g cm~^ (black) and 

= 2 g cm“^ (blue). 

level (which confirms the ®®Ni mass, the density, and the abundance 
stratification) and in its development over time (which is a result of 
photon diffusion and again tests density, ®®Ni mass, and abundance 
distribution). The agreement in the rising part of the light curve is 
remarkable. This would have been worse had we adopted a shorter 
rise time. More mixing out of ®®Ni might have been required in that 
case, bu t this would be in co nflict with results from the early-time 
spectra dMazzali et aLlbOl^ . The main shortcomings of our syn¬ 
thetic light curve are the slight underestimate of the flux just after 
maximum brightness (see Figure [8] inset) and the overestimate at 
100-200 days. The discrepancy just after peak may indicate that we 
are underestimating the opacity in the inner ejecta, which is quite 
possible given our rather coarse grey treatment. The discrepancy 
at day 100-200 may be due to the opposite problem, but also to a 
less well determined correction for flux outside the optical bands, 
which may shift the bolometric light curve of SN 201 Ife upward. 


5 DISCUSSION 

Our results appear to confirm that the density structure derived by 
IMazzali et al.l ( 12014l) can appropriately describe SN201 Ife when it 
is extended to low velocities under the assumption that the explo¬ 
sion ejected approximately one Chandrasekhar mass of material. 
They also strongly favour the presence of a significant amount of 
®^Fe, as long as it is located at the lowest velocities. The mass of 
®®Ni produced by SN201 Ife is ~ 0.45 M©, in agreement with es¬ 
timates from the light curve, if ®®Ni is located above ®'*Fe. On the 


other hand, if ®®Ni is located centrally, the mass which is necessary 
to power the late-time spectra is significantly smaller. The amount 
of stable ®®Ni required to match the observed emission near 7400 A 
is very small, ~ 0.06 Mq, because that emission feature is domi¬ 
nated by [Fe It] and [Ca It] lines. 

While it is difficult to determine the mass with high precision 
(even in the case of SN 201 Ife the uncertainties in luminosity, com¬ 
bined with the intrinsic approximations of the model, mean that we 
cannot claim better than 20% accuracy), we note that the kinetic 
energy of model p-llfe matches very well the expected Ek yield 
based on the nucleosynthesis. Applying an equation relating the 
yield of different species to the Ek produced in a Chandrasekhar- 
mass SN la explosion, 

Ek = [1.56 M(®®Ni) + 1.74 M(stableNSE) + 

1.24 M(IME) - 0.46] X 10®^ erg, (1) 

dWooslev et al.ll200^ . where the constant term is the binding en¬ 
ergy of the WD, and inserting M(®®Ni) = 0.47 M©, M(stable 
NSE) = 0.24 Mq, and M(IME) = 0.41 M©, we obtain Ek = 
1.20 X 10®^ erg, in excellent agreement with the energy of p-1 Ife, 
1.32 X 10®^ erg. This suggests that a weak delayed detonation 
of a near-Chandrasekhar-mass white dwarf is a consistent model 
for SN 201 Ife. The fact that the emission-line width of SN 201 Ife 
is low for its light-curve shape (the width of the strongest lines 
is very similar to those of SNe 1992A and 2004eo, which had 
Ami 5 (i?) « 1.4 mag) also indicates that the specific energy of 
SN201 Ife was somewhat lower than normal. 

Other points are worthy of notice, as follows. 

(1) 1D-3D. Our ID model p-llfe resembles delayed- 
detonation models th at were computed in ID, such as those of 
llwamoto et al.l ( Il999l) . ID delayed-detonation models predict that 
n-rich NSE species are located deeper than ®®Ni. Recently, a 3D 
delayed-detonation model has bee n published which sugg ests that 
®®Ni is located below stable Ee dSeitenzahl et alj |201^ . A test 
with a simplistic reproduction of that situation does not lead to an 
equally good synthetic spectrum. There seems to be no reason why 
our ID synthetic spectra should be biased in favour of ID explosion 
models, as these are computed completely independently. 

(2) Mass of progenitor. An accurate estimate of the mass does 
not seem to be possible to better than ~ 20%, because of uncertain¬ 
ties related to data calibration, distance and reddening, as well as to 
the modelling process. However, our models require the presence of 
n-rich NSE species at the deepest velocities in order to achieve an 
ionisation regime that allows the ratio of the strongest [Ee III] and 
[Ee II] lines to be reproduced. This is natural in SD models where 
the progenitor WD is approaching the Chandrasekhar mass and 
has a high central density. It is not in edge-lit sub-Chandrasekhar- 
mass models, or in merger models, where in both cases the cen¬ 
tral density is lower. These two classes of events may give rise 
to rather distinct spectroscopic signatures in the nebular phase , 
which could explain some s ubtypes of SNe la dMazzali et aljTOlU : 
IMazzali & HachingejTOlj) . On the other hand, we cannot exclude 
a mass somewhat different from the Chandrasekhar mass, or a dif¬ 
ferent density profile. In particular, the nebular lines are narrow 
for a SN with the light-curve shape of SN201 Ife {Amis{B) « 
1.1 mag). The model we have presented here seems to show a high 
degree of consistency. It needs to be tested whether first-principle 
3D models can explain SN 201 Ife. T his could be done using codes 
such as those of iMaeda et al.l (l2006l) . 

(3) The Ni line. If the line near 7400 A is indeed [Ni II] 
AA7380, 7410, a small amount of ®®Ni is sufficient to reproduce it. 
A possible blueshift (~ 1100 km s~^: iMcClelland et alj|2013l) is 
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observed with respect to the rest-frame wavelength of the lines, af¬ 
ter correcti on for cosmological redshift. In the aspherical explosion 
scenario of iMaeda et alJ ll 2010 h . this is consistent with SN 201 Ife 
being a low velocity gr adient (LVG) SN (t) = 35 ± 5 km s“^d“^; 
[McClelland et al.ll2013h . It should be noted that the [Ni II] line is 
blended with the emission near 7300 A, which is dominated by 
[Ca II] and [Fe II] lines. Some of the blueshift of the Ni line may 
simply be the result of this blending. 

One of the ingredients of the 3D model sketched by 
[MaedaetalJ ( l 2010 l) is that the centrally located n-rich NSE species 
are not excited by 7 -rays, and that this would be revealed by flat- 
topped profiles of some emission lines. In particular, in the case of 
SN 2003hv, evidence of flat tops was foun d in the NIR spectrum 
jMotohara et alj|200^ : iCerardv et alj[2007ll . However, our models 
show that it is not difficult for inwardly propagating 7 -rays emit¬ 
ted from a ®®Ni zone that surrounds a nonradioactive central core 
to excite material located at the lowest velocities. This is indi¬ 
cated by the sharp peak of the [Ni II] emission. Also, in the case 
of SN 201 Ife, the NIR spectrum does not seem to show any flat- 
topped profiles. This suggests that not only is ®®Ni irradiated, but 
®^Fe is as well. In particular, the line at 1.26 /rm is dominated by 
a single line, [Fe II] 1.257 fim, and hence is the perfect candidate 
for a flat top, but we see no evidence for one. Incidentally, a flat 
top in that line d oes not seem to be v isible in the spectrum of 
SN 2003hv either dMotohara et al.|[200^ . an d the emission coul d 
be reproduced quite well with a ID model dMazzali et aklbOlTh . 
The other line for which a flat top was observed in SN 2003hv is at 
1.6 fim. This line is a blend of [Fe II] and [Si II]. It is not impos¬ 
sible that the observed flat top is the result of the blending of these 
lines and the relatively low signal-to-noise ratio of the NIR spec¬ 
trum of SN 2003hv. In SN 201 Ife no flat top is present in this fea¬ 
ture, but SN 201 Ife was more luminous than SN 200 3hv, and there¬ 
fore i t probably produced less Si and other IMEs dMazzali et alj 
l2007h . So, if the flat tops in SN 2003hv are real, the density con¬ 
trast between the n-rich core and the surrounding ®®Ni shell was 
perhaps very large, so that 7 -rays could not penetrate. This may 
be unlikely if SN 2003hv wa s a sub-Chandrasekhar-mass explosion 
dMazzali & Hachingei[|2012h . but it may be worth testing. 

(4) Hydrogen? N o narrow Ha emissi on is seen in the neb¬ 
ular spectra (see also iShappee et al.l [201 3h . An emission line is 
observed at a similar wavelength, and the model indicates that it 
may be [Ee II] A6559, but the synthetic line does not match the ob¬ 
served one in flux. Several other weak lines are poorly reproduced. 
They are mostly attributed to Ee lines because of wavelength co¬ 
incidence (see Fig.|5}. Their incorrect strengths may be caused by 
poorly known collision strengths. However, if the distribution of 
H stripped fr om a companion was not as narrow as p redicted by 
some models dMarietta et al.ll2Q0^ : IPakmor et al.ll201 ih . Ha could 
make a contribution to the observed emission. The observed line 
has similar width as all other lines. 


6 CONCLUSIONS 

We have presented ID models for a series of optical and one NIR 
spectrum of the nearby, spectroscopically normal SN la 201 He. We 
have shown that if we use the density dis tribution derived fro m the 
early-time optical/UV spectra (“g-1 Ife’k lMazzali et alj|2014ll . and 
extend it to low velocities assuming that the ejected mass equals 
the Chandrasekhar mass, we can obtain a satisfactory match of 
the main features of the late-time spectra using a ®®Ni mass of 
~ 0.45 M© and a stable Fe mass of ~ 0.23 M©. Stable Fe should 


be mostly located in the central, lowest-velocity region of the ejecta 
for a best match to be ob tained. This is consis tent with ID delayed- 
detonation models (e.g.. ilw amoto et alJl999ll . but not with a recent 
3D model dSeitenzahl et al.l2013L but see Gamezo et al. 2005). The 
derived density and abundance distributions were used to compute 
a bolometric light curve, in the spirit of abundance tomography. 
The result compares favourably with th e pseudo-bolometric light 
curve of SN20life ( jPereira et al. I I2OI3I) . lending further credence 
to the model proposed here. 

A significantly sub-Chandrasekhar model seems to be ruled 
out because the low central density and the ensuing lack of newly 
synthesised ®^Fe lead to a high ionisation and incorrect ratios of 
the strongest [Fe II] and [Fe III] lines. Merger models that involve 
WDs significantly below the Chandrasekhar mass may have simi¬ 
lar characteristics, and would then be equally disfavoured on these 
grounds. 

A weak line near 7400 A can be reproduced as [Ni II], requir¬ 
ing the presence of a small amount of ®®Ni (~ 0.06 Mq). The line 
shows a blueshift of ~ 1100 km s~^ with respect to its rest wave¬ 
length, in agreement with the LVG nature of SN 201 He. The central 
Fe region is effectively excited by 7 -rays from the ®®Ni region that 
surrounds it, preventing the formation of flat-top emission profiles 
even for rather isolated Fe lines. 

The evidence concernin g the nature of SN 201 He remains am¬ 
biguous teopke et alj|2012h . Spectral fits seem to favour a model 
resem bling a ID delayed detonation, but with lower energy (this 
paper; [Mazzali et al .I2OI4I) . On the other hand, much observational 
evid ence, including the lack of hydrogen in the late-time spec- 
tra liShap pee et al .] [2013[) . the lack of X-ray or radio detections 
jHoresh e t al.l20i ai. the a bsence of early-time extra flux in the light 
curve Isi oom et a l][ 2012 h . and the lack of an imaged companion 
jLi etal.[[201lh . does not support SD models with a main sequence 
or red giant companion. On the other hand, models that do not reach 
high central densities and do not produce stable NSE isotopes, in¬ 
cluding significantly sub-Chandrasekhar models (such as might re¬ 
sult from the accretion of helium from a nondegenerate He compan¬ 
ion with small radius) or mergers of WDs whose individual mass is 
well below the Chandrasekhar mass, are not favoured by late-time 
models. Additionally, no evidence for any residual He was se en in 
the early-time spectra l lParrent et al .120 13 ; [Mazzali et alj2014h . On 
the other hand, a marginally sub-Chandrasekhar model may be able 
to reproduce both the low energy and the presence of stable NSE 
isotopes. Alternatively, a marginally super-Chandrasekhar model 
may also display all observed features. 

In summary, spectroscopic evidence seems to favour for 
SN 201 He a model with a high central density, matching the struc¬ 
ture of a low-energy delayed detonation of a Chandrasekhar-mass 
WD in a SD system. The model we used has Ek = 1.20 x 10®^ erg. 
Currently, ID models offer a better match to the data than 3D mod¬ 
els, but this should be taken with some caution because ID mod¬ 
els are intrinsically simpler. Further work on this and other models 
would be useful to bring them into agreement with all observed 
properties of SNe la. 
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